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Abstract 



The correct understanding of the nature and dynamics of interatomic magnetic 
interactions in sohds is fundamentally important. In addition to that it allows to 
address and solve many practical questions such as stability of equilibrium mag- 
netic structures, designing of magnetic phase diagrams, the low-tempcraturc spin 
dynamics, etc. The magnetic transition temperature is also related with the be- 
havior of interatomic magnetic interactions. One of the most interesting classes of 
magnetic compounds, which exhibits the rich variety of the above-mentioned prop- 
erties in the transition-metal oxides. There is no doubts that all these properties 
are related with details of the electronic structure. In the spin-density-functional 
theory (SDFT), underlying many modern first-principles electronic structure meth- 
ods, there is a certain number of fundamental theorems, which in principles provides 
a solid theoretical basis for the analysis of the interatomic magnetic interactions. 
One of them is the magnetic force theorem, which connects the total energy change 
with the change of single-particle energies obtained from solution of the Kohn- 
Sham equations for the ground state. The basic problem is that in practical im- 
plementations SDFT is always supplemented with additional approximations, such 
as local-spin-dcnsity approximation (LSDA), LSDA -|- Hubbard U, etc., which are 
not always adequate for the transition-metal oxides. Therefore, there is not perfect 
methods, and the electronic structure we typically have to deal with is always ap- 
proximate. The main purpose of this article, is to show how this, sometimes very 
limited information about the electronic structure extracted from the conventional 
calculations can be used for the solution of several practical questions, accumulated 
in the field of magnetism of the transition-metal oxides. This point will be illus- 
trated for colossal-magnetoresistive manganites, double perovskites, and magnetic 
pyrochlores. We will review both successes and traps existing in the first-principle 
electronic structure calculations, and make connections with the models which cap- 
ture the basic physics of the considered compounds. Particularly, we will show what 
kind of problems can be solved by adding the Hubbard U term on the top of the 
LSDA description. It is by no means a panacea from all existing problems of LSDA, 
and one should clearly distinguish the cases when U is indeed indispensable, play a 
minor role, or may even lead to the systematic error. 
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1 Introduction 

The first-principles electronic structure calculations play a very impor- 
tant role in the exploration of magnetism. They have been successfully 
applied for various types of metallic compounds. The transition-metal ox- 
ides (TMO), however, take a very special place in this classification and 
typically regarded as a counter-example where the first-principles calcu- 
lations either experience serious difficulties or simply do not work. Such 
an extreme point of view has of course a very serious background because 
most of modern computational techniques are based on the spin-density- 
functional theory (SDFT), designed for the ground state, which is typically 
supplemented with the local-spin density approximation (LSDA) for the 
exchange-correlation interactions. The latter is based on the homogeneous 
electron gas theory, and therefore is very different from the localized-orbital 
hmit, which was originally adopted for the description of TMO [U [21 El HI • 
Therefore, traditionally there was a big gap in the understanding of in- 
teratomic magnetic interactions in TMO basing on these model arguments 
and the first-principles electronic structure calculations jSl E] , which ap- 
peared only in the middle of 1980s and were regarded as very challenging 
at that time. 

Since then, two different standpoints have certainly became closer. It 
is true that due to complexity of the problem of exchange and correlations, 
even now there is no perfect computational methods for the transition- 
metal oxides. However, it also becomes increasingly clear that the model 
analysis of the problem should use results of first-principles electronic struc- 
ture calculations, at least as a starting point. In many cases, a puzzling 
behavior attributed to the fanciful correlation effects can be naturally ex- 
plained by details of the reahstic electronic structure. 

The purpose of this article is to make a link between general formula- 
tion of SDFT and models of interatomic magnetic interactions for TMO. 
We will try to show that even in the present form the electronic struc- 
ture calculations can play a very important role in the understanding of 
magnetic properties of various TMO, despite many limitations inherent to 
LSDA and some of its refinements. 

In Sec. El we will summarize main results of SDFT and formulate the 
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magnetic force theorem which presents the basis for the analysis of mag- 
netic interactions in solids. The connection of this theorem with some 
canonical models for interatomic magnetic interactions will be illustrated 
in Sec. El In Sec. HI we will present a critical analysis of interatomic mag- 
netic interactions in MnO basing on several available methods of electronic 
structure calculations. In Sec. El we will consider some practical problems 
related with the colossal-magnetoresistive (CMR) manganites, double per- 
ovskites Sr2FeM06 {M= Mo and Re), and magnetic pyrochlores 742M02O7 
{A= Y, Gd, and Nd). 

The insulating character of many TMO presents one of the most in- 
teresting and controversial problems. Particularly, it is well known that 
LSD A frequently underestimates or even fails to reproduce the energy gap, 
which is formally the excited state property. Is it possible that even in 
this case, it can provide a physically meaningful description for the param- 
eters of the ground state, such as the interatomic magnetic interactions? 
In this context, we will consider the role played by the on-site Coulomb 
interaction U on the top of the LSDA electronic structure and argue that 
one should clearly distinguish the behavior of band insulators, which is 
governed by the double exchange (DE) physics j2l EI , and in principle can 
be accounted for by itinerant models, in the spirit of LSDA, from the be- 
havior of Mott insulators, where the Coulomb U is indeed indispensable 
in order to suppress (in this case) spurious DE interactions and unveil the 
completely different physical behavior governed by the superexchange (SE) 
interactions |3l HI IS|- A brief summary and perspectives will be outlined in 
Sec. El 

2 Theory of magnetic interactions in solids 
2.1 Spin-Density- Functional Theory 

The modern way to approach the problem of interatomic magnetic inter- 
actions in solids is based on the SDFT [9j , which states that the magnetic 
ground state of the A^-electron system can be obtained by minimizing the 
Hohenberg-Kohn total energy functional 

E[m] =ro[m] + ^xc[m] (1) 
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(Tofm] being the kinetic energy of a non-interaction electron system and 
E'xcfm] being the exchange-correlation energy) with respect to the spin- 
magnetization density m(r).^ 

The direct implementation of the SDFT is hampered by the fact that 
the functional dependencies To[m] and £'xc[m] are generally unknown. In 
the case of Tofm], the problem is resolved by introducing an auxiliary 
system of one-electron orbitals {ipi{r)}, and requesting the kinetic energy 
(in Rydberg units), 



N 

z 

i=l 

the spin-magnetization density, 



^oH = E /«5rV^J(r) (-V^) ij,{r] 



m(r) = E^l(r)^T^^(r) (2) 

(<T being the vector of Pauli matrices), and the total energy to coincide 
with the same parameters of the real many-electron system in the ground 
state. Then, the minimization of £"[111] with respect to m(r) is equivalent 
to the self-consistent solution of single-particle Kohn-Sham (KS) equations 
for {V^z(r)}: 

-V^ + fT-B{r)]ilj,{r)=eMr), (3) 
where the effective magnetic field is given by 

B(r) = -^£;xc[m]- 

The exchange-correlation energy functional, £'xc[ni], is typically taken 
in an approximate form. Several possible approximations along this line 
are listed below. 

• The local-spin-density approximation. In this case, the explicit de- 
pendence of E'xc on m(r) and the electron density n(r), 

Exc[n,in] = J drn{r)£xc [n{r), |m(r)|] , (4) 



^For the sake of simplicity, we have dropped in Eq. Q the electron density, n{r), and all terms which 
depend only on n{r), though this dependence is implied as well as the minimization of E with respect to 
n{r). 
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is borrowed from the theory of homogeneous electron gas. Concep- 
tually, LSDA is similar to the Stoner theory of band magnetism [lOj. 
Due to the rotational invariance of the homogeneous electron gas, Exc 
depends only on the absolute value of m(r). 

• The local-density approximation plus Hubbard U (LDA+?7) ap- 
proach pn [121 US]- This is a semi-empirical approach, the main idea 
of which is to cure some shortcomings of the LSDA description for the 
localized electron states by replacing corresponding part of E'xc[^, m] 
in LSDA by the energy of on-site Coulomb interactions, in an analogy 
with the multi-orbital Hubbard model. The latter is typically taken 
in the mean-field Hartree-Fock form: 

-E'xc[^ ] = X! (^71737274 ~ ^71737472)^7172^7374' 

{7} 

where ^^71 737274 =(71 73 1 — 1 7274) are the matrix elements of the on-site 
Coulomb interactions, which are assumed to be renormalized from the 
bare atomic values by interactions with other (itinerant) electrons and 
by correlation effects in solids. The electron and spin-magnetization 
densities for the localized states are represented by corresponding el- 
ements of the density matrix n'^=||n!^^^2 1| in the basis of atomic- like 
orbitals {7} at the site r. Because of this construction, the LDA+[/ 
approach is basis-dependent and typically implemented in the linear 
mufhn-tin orbital method jl4j . 

• The optimized effective potential (OEP) method jlSlIIB] is an attempt 
of exact numerical solution of the KS problem, which does not rely 
on the local-spin-density approximation for E'xcf^jin]- In this case, 
the eigenfunctions {?/^i(r)} and eigenvalues {si} obtained from the 
KS equations (01) with some trial potential are used as an input for 
total energy calculations beyond the homogeneous electron gas limit. ^ 
The parameters of such potential are requested to minimize the total 
energy. 

^Of course, practical implementations of the OEP scheme require some approximations for the 
exchange-correlation energy. Typically, it is either Hartree-Fock or the the random-phase approxima- 
tion, underlying the GW method [T7)'). 
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2.2 Magnetic Force Theorem 



The basic idea behind analysis of interatomic magnetic interactions in 
sohds is to evaluate the total energy change AE=E[RQeGs]—E[eGs] caused 
by small non-uniform rotations of the spin-magnetization density near the 
ground state, where eGs(r)=niGs(i")/|niGs(r)| is the direction of the spin- 
magnetization in the ground state and Re is the three-dimensional rotation 
by the (small) angle ^(r), which depends on the position r in the real space: 

^0(r)eGs(r) = eGs(r) + [6>(r) x eGs(r)] - ^6>2(r)eGs(r). 

Thus, characterizes the local stability of the magnetic ground state with 
respect to non-uniform rotations of the spin-magnetization density. 

Very generally, the problem can be solved by minimizing the con- 
strained total energy functional: 

Eh[m] = ToH + ExcH - / drh{r) ■ [e(r) - Rg^^^eGsir)] , 

where the constraining field h(r) plays the role of Lagrange multipliers and 
enforces the requested distribution of the spin-magnetization density in the 
real space. 

The magnetic force theorem states in this respect: in the second order 
of 6{r) the total energy change is solely determined by the change of the 
KS single-particle energies 



N 

AE=Y: {e, [ReB, ReuiGs] - [B, niGs]) + O {O^) . (5) 

The eigenvalues Si RoB^Rquigs , corresponding to the rotated effective 
field RqB and the ground state spin-magnetization density i^^niGs, can be 
expressed in terms of expectation values of the KS Hamiltonian, 

£, [ReB, ReuiGs] = J dr4[ReinGs] (-V^ + o- • ReB) V^.fi^eniGs], 



with {t/jilRouiGs]} yielding i^^niGs after substitution into Eq. 

The theorem can be reformulated in a different way: the effect of the 
longitudinal change of m(r) on Tofm] and £'xc[ni] caused by the self- 
consistent solution of the KS equations (01) with the external magnetic 
field h(r) is cancelled out in the second order of 0(r). 



7 



The theorem can be proven rather generally provided that the 
exchange-correlation energy functional obeys the following condition jlB] : 

E^cH = Ey^ciReT^]- (6) 

In LSDA, it immediately follows from Eq. (jlj). More generally, Eq. (Ej) can 
be regarded as the fundamental gauge-symmetry constraint, which should 
be superimposed on the admissible form of the exchange-correlation energy 
functionals fm|.^ 

The magnetic force theorem has very important consequences: 

• Generally, the KS eigenvalues {si} have no physical meaning and can- 
not be compared with the true single-particle excitations (that is typ- 
ically the case for many spectroscopic applications). In this respect, 
the magnetic excitations present a pleasant exception, thanks to the 
magnetic force theorem. 

• In principle, the knowledge of the effective KS potential alone is suf- 
ficient to calculate the total energy difference. It allows to get rid of 
heavy total energy calculations (particularly, for the OEP method) 
without any loss of the accuracy. 



2.3 Practical Implementations of the Magnetic Force Theorem 

In most cases, practical applications of the magnetic force theorem are 
based on relaxed constraint conditions in comparison with the ones given 
by Eq. (Ej). Namely, the second condition requesting the KS eigenvalues to 
correspond the rotated spin magnetization density is typically dropped and 
ReB^RouiQs is replaced by £i RqQ obtained from the KS equations 
with the rotated field i^^B. This considerably facilitates the calcula- 
tions, though at the cost of a systematic error, which was recently discussed 
by Bruno [20j.^ We will return to this problem in Sec. I5.1.il where will give 

•^Since m(r) is given by Eq. JJjl, the three-dimensional rotation m(r)^i?e(r)m(r) is equivalent to 

the unitary transformation of the KS orbitals tpi{r)^Us [^(r)] i^iii^), where Us [6{r)]=cxp [^cr ■ 0(r)] is 
the 2x2 rotation matrix in the spin subspace. Then, the property ^ can be easily proven for many 
other functionals, for example the ones based on the Hartree-Fock approximation and underlying the 
rotationally invariant LDA+C/ and OEP ^H] methods. 

* Since B (r ) = -g^^^ Exc [m] and Exc [m] satisfies Eq. , the rotation of the spin- magnetization density 

mGs(r)^i?e(r)HiGs(r) results in similar rotation of the effective field B(r)^i?e(r)B(r). Therefore, at the 
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some estimates of this error and discuss some physical implications which 
may be related with the additional rotation of the spin-magnetization den- 
sity for the transition-metal oxides. In this section we will review some 
practical schemes of calculations, which ignore these effects. Basically, 
they are two. 

The first one is based on the perturbation theory expansion for the 
Green function ^Tj 



and its projections 



£ — £i + i6 



on the majority (f ) and minority (j) spin states (Tr5' being the trace over 
the spin indices). In the second order of 0{r), the total energy change can 
be mapped onto the Heisenberg model 



AE = -^1 dr I dr'J{r,r') [i?0(r)eGs(r) • Ro{r')eGs{r') - eGs(r) • eGs(rO] 



The parameters of this model are given by 



(7) 



J(r, r') = -Im f'^ deG^r, r', £)B{r')G''{r' , r, £)B{r), (8) 

where £f is the Fermi energy corresponding to the highest occupied KS or- 
bital. All practical calculations along this line are typically performed on a 

input of the KS equations ((SJ the direction of the spin-magnetization is consistent with the direction of 
the effective field. However, the direction of the new magnetization obtained after the first iteration can 
be canted off the initial distribution prescribed by the matrix Rg. This is precisely the source of the error, 
and strictly speaking the magnetic force theorem cannot be proven for this relaxed constrained condition. 
In order to correct this error, Bruno [20] explicitly considered the constraining field h(r), which can be 
estimated using the response-type arguments. This field does not explicitly contribute to the expression 
© for the total energy change. However, it does modify the KS orbitals {ipi}, which should be used 

in order to evaluate Si i?eB, .RemGs • Although the analysis undertaken by Bruno is certainly valid 

and the finding is important, we believe that there is some confusion with the terminology, because the 
magnetic force theorem itself is correct, as long as it is formulated in the form of Eq. ||SJ) The main 
criticism by Bruno j20j is devoted to the practical implementations of this theorem. 

^Here, the mapping onto the Heisenberg model is a general property, whereas the from of Eq. ijSJ 



corresponds to the additional approximation Si 



Si i?6(B for the KS eigenvalues. 
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discrete lattice, assuming that all space is divided into atomic regions spec- 
ified by the site indices {r}, so that the angles 6{r) depend only on r, and 
neglecting the effects caused by rotations of the spin-magnetization density 
on the intra-atomic scale. ^ In the discrete version, G'^'-'-(r, r', e) is expanded 
in the basis of atomic-like orbitals {7} centered at different atomic sites 
(for example, in the nearly-orthogonal LMTO representation |[I3j)- The 
integral over atomic regions around r and r', Jrr'^SrdY Sr' dr'J{r, r'), de- 
fines conventional parameters of interatomic magnetic interactions in the 
real space. 

The main idea of the second approach is to calculate the energy of the 
collective spin excitation corresponding to the frozen spin wave with the 
vector q. The method, which is called the frozen (or adiabatic) spin-wave 
approximation is designed for the discrete lattice [2S|- In this context, the 
adiabaticity means that the directions of the spin-magnetization can be 
regarded as the "slow" variables, so that for each configuration specified 
by Rg the "fast" electronic degrees of freedom have enough time to fol- 
low this directional distribution of the spin-magnetization density. The 
rotation matrix Rq is requested to transform the coUinear distribution of 
the spin magnetic moments in the ground state, mT-=/T-<ir mGs(r), to the 
spin spiral: mT-^(cosq • Tsin6'T-,sinq • Tsin^T-,cos^T-)|mT-|. If mT-||z, the 
rotation angles are given by 

^T- = (— sin q • T, cos q • r, 0) 9t 

{6r being the cone-angle of the spin wave). The KS equations for the spin- 
spiral configuration of the effective field RqB can be solved by employing 
the generalized Bloch transformation [24j. This gives the total energy 
change AE'(q, 0) corresponding to the "excited" spin-spiral configuration 
with the vector q. In the the second order of 0, this energy change can 
be mapped onto the Heisenberg model. The mapping provides parameters 
of magnetic interactions in the reciprocal space, Jq, which can be Fourier 
transformed to the real space. 

^The non-coUinear distribution of the spin-magnetization density on the intra-atomic scale is an in- 
teresting and so far very imperfectly understood phenomenon |22| . 
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3 Relation with the Models of Magnetic Interactions 



The application of Eq. (IHI) goes far beyond the standard electronic struc- 
ture calculations. It is rather general, and many canonical expressions for 
magnetic interactions in solids can be derived starting with Eq. (jHI). Here 
we would like to illustrate this idea by considering two model examples. 



3.1 Double Exchange and Superexchange Interactions in Half- 
Metallic Ferromagnetic State 

Consider the ferromagnetic (FM) chain of atoms, described by the Hamil- 
tonian 

which is an analog of the KS Hamiltonian (01) on the discrete lattice, where 
matrix elements of the kinetic energy Iq (the transfer integrals) are re- 
stricted by the nearest neighbors, and B is the effective field polarizing 
the conduction electrons parallel to the z axis. The half-metallic behav- 
ior implies that to<B and the electron density n<l (which corresponds 
to the partial population of the j-spin band in Fig. [T]). The problem can 



o 



Q 




B-t B B+t„ 



-B-t„ -B -B+t 



h 



Figure 1: Schematic density of states for the half-metalhc ferromagnetic chain. 

be easily solved analytically, and by expanding gI.^,{£) in Eq. (jHI) up to 
the second order ofto/B one can obtain the following contributions to the 
nearest-neighbor (NN) exchange coupling 



= — sm Trn 
2n 



to 



(9) 
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and 

= --^(7m + -sin27rn). (10) 

is the FM double exchange interaction, which is proportional to to. 
is the antiferromagnetic (AFM) superexchange interaction. For the half- 
filled band at n=l the system is insulating. Then, J^=0, while Eq. ( jTUI) 
leads to the standard expression for the SE interaction at the half-filling: 

More generally, and can be expressed through the moments of 
local density of states, and is the measure of the kinetic energy of the 
fully spin-polarized half-metalhc states jSHI • Originally, the concept of the 
double exchange was introduced for the metallic phase of CMR mangan- 
ites [7]. However, the phenomenon appears to be more generic and one of 
the most interesting recent suggestions was that the same mechanism can 
operate in the insulating state, provided that the system is a band insula- 
tor \2§\ ■ This substantially modifies the original view on the problem and 
combine two seemingly orthogonal concepts, one of which is the DE physics 
and the other one is the insulating behavior of some CMR manganites. We 
will return to this problem in Sec. 15. IL 



3.2 Superexchange Interaction via Oxygen States 

Consider the interaction between two magnetic (transition-metal) site, 
which is mediated by the non-magnetic oxygen states (Fig. [21 the situ- 
ation is rather common for the insulating transition-metal oxides). It is 
assumed that the sphtting between the j- and j-spin states of the tran- 
sition metal sites is 25. A describes the relative position of the oxygen 
states relative to the transition-metal states (the so-called charge-transfer 
energy). It is further assumed that the occupied states are Ti(f), T2(i), 
and 0(1,1), and the unoccupied states are Ti{i) and T2(t). 2B and A are 
the largest parameters in the problem, so that the transfer interactions t 
between the transition-metal and oxygen sites can be treated as a pertur- 
bation starting with the atomic limit. Corresponding matrix elements of 
the Green function connecting the sites Ti and T2 are given by 

GUe) = gUs) = t\e + B)-\e -B + A)-\e - B)-\ (11) 
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Figure 2: Positions of atomic levels illustrating the superexchange interaction between 
transition-metal sites mediated by the oxygen states. 

Then, using Eq. (IHI) one can obtain the well-known expression for the SE 
interactions mediated by the oxygen states f27\-J 



which describes the shift of the poles of the Green function ( fTTI) of the 
occupied states, located at £=—B and £=B—A, due to the interaction 
with the unoccupied states located at £=B. 

Considered examples of DE and SE interactions show that Eq. (jSI) is 
rather universal and can be regarded as the starting point for the general 
analysis of interatomic magnetic interactions in solids. The applications 
can be very wide and cover, for example, the theory of RKKY interac- 
tions [2B|, the effects of the interatomic Coulomb interactions on the SE 
coupling [29j, etc. 

^Here we have considered the AFM configuration of the sites Ti and T2 (Fig. The assumption 
is not important and absolutely the same expression for can be obtained by starting with the FM 
configuration. This means that in the case of superexchange, the mapping of the total energy change onto 
the Heisenberg model is universal. This is not a general rule, and other types of magnetic interactions do 
depend on the state in which they are calculated. This dependence is related with the change of electronic 
structure which enters the expression for the magnetic interactions © through the matrix elements of 
the Green function. As we will see in Sec. 15.11 the form of the DE interactions strongly depends on the 
magnetic state and this dependence plays a crucial role in the physics of CMR manganites. 
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4 Electronic Structure of MnO from the viewpoint 
of Interatomic Magnetic Interactions 

The manganese monoxide provides an excellent opportunity to test avail- 
able methods of electronic structure calculations for the transition-metal 
oxides. While numerous spectroscopic techniques deal mainly with the 
excited state properties, which cannot be compared directly with results 
of electronic structure calculations designed for the ground state, the in- 
teratomic magnetic interactions are the ground-state properties and, in 
principle, should be accounted for by these calculations. 

The interatomic magnetic interactions in MnO are well studied exper- 
imentally. The magnetic behavior of MnO can be described by the simple 
Heisenberg model including first (Ji=— 4.8 meV) and second (J2=— 5.6 
meV) neighbor interactions |30j.^ Since MnO is a wide-gap insulator, it 
is clear that both interactions originate from the superexchange mecha- 
nism p[ (S^ . 

One can easily design the proper model for the electronic structure of 
MnO using the following arguments [TB|. Since the electronic configura- 
tion of the Mn atoms in MnO is close to 3d^3dP^, both from the viewpoint 
of the model valence arguments and the electronic structure calculations, 
the distribution of the spin-magnetization density near Mn sites is nearly 
spherical. Therefore, the only parameter of the effective magnetic field at 
the Mn sites we need to worry about is 5, which controls the splitting 
between occupied and unoccupied 3(i| states. Another important pa- 
rameter of the electronic structure is the charge-transfer energy A, which 
controls the relative position of the 3(i| and the oxygen 2p states. Then, 
using an analogy with the model for SE interactions considered in Sec. 13. 2L 
one can try to parameterize the effective KS potential in terms of B and 
A (basically in the spirit of LDA+[/ for the half-filled 3d shell), and to 
find these two parameters by fitting two experimental parameters of in- 
teratomic magnetic interactions, Ji and J2, provided that the kinetic and 
non-magnetic parts of the effective KS Hamiltonian are well described on 

^According to the definition ITJ of the Heisenberg model, the experimental parameters have been 
multiphed by 5^ = (5/2)^. We do not consider the trigonal distortion of the cubic lattice caused by the 
exchange striction effects 
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the level of LDA. 

The fitting, which is very straightforward in the nearly-orthogonal 
LMTO basis [HI, [TH], yields the following parameters: i?~5.3 eV and 
A~10.7 eV.^ The corresponding density of states is shown in Fig. El 
together with results of the LSDA, LDA+t/, and OEP calculations.^^ 
Corresponding parameters of magnetic interactions are listed in Table [H 
Unfortunately, none of the existing techniques can treat the problem of in- 
Table 1: Parameters of magnetic interactions in MnO obtained in LSDA, OEP and 
LDA+f/ in comparison with experimental data |3Uj . 



teratomic magnetic interactions in MnO properly. Nevertheless, the anal- 
ysis of the electronic structures shown in Fig. El allows us to elucidate the 
basic problems of the existing methods which lead to substantial overesti- 
mation of Ji and J2 [THj . 

• LSDA underestimates both B and A. In this sense, the overestima- 
tion of interatomic magnetic interactions is directly related with the 
underestimation of the band gap. However, the occupied part of the 
spectrum is reproduced reasonably well by the LSDA calculations. 

• LDA+[/ underestimates A. This is a very serious problem inherent 
not only to LDA+t/ but also to some subsequent methodological de- 
velopments for the strongly correlated systems like LDA+DMFT ap- 
proach [211 , which provides a solid basis for the description of Coulomb 
interactions at the transition-metal sites, but ignore other parameters 
of electronic structure, which may play an important role for TMO. 

^The on-site Coulomb repulsion U for the 2>d electrons can be estimated from 2B by using the connec- 
tion 2B=^{U + AJ)iiMn which holds for the mean-field Hartree-Fock solution of the degenerate Hubbard 
model (/iMn— 4.84 hb being the local magnetic moment). Then, using J~0.8 eV for the intra-atomic 
exchange coupling, which is not sensitive to the crystal environment in solid one obtains J7~8 eV. 

^"AU OEP calculations discussed in this section have been performed using the one-electron effective 
potentials obtained by T. Kotani 



Ji (meV) J2 (meV) 



LSDA 
OEP 



-13.2 -23.5 

-8.9 -12.0 

-5.0 -13.2 

-4.8 -5.6 



LDA+f/ 
Expt. 
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Figure 3: Electronic structure of MnO, which fits the experimental parameters of mag- 
netic interactions together with results of first-principles calculations in LSDA, OEP, and 
LDA+C/. 

• The relative position of the 3d and 2p states is well reproduced by 
the OEP method. However, the width of the unoccupied 3d band is 
strongly overestimated. This might be related with the local form of 
the effective one-electron potential in the OEP approach. As the re- 
sult, the parameters of interatomic magnetic interactions are strongly 
overestimated. 

Thus, the quantitative description of interatomic magnetic interactions 
in TMO by the first-principle electronic structure calculations is largely un- 
resolved and challenging problem. The analysis of MnO suggests that one 
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possible direction along this line could be to borrow the form of the one- 
electron potential from the LDA+t/ approach and try to optimize param- 
eters of this potential by using the ideas of variational OEP method [.18j . 
In any case, the development of the first-principle techniques, which could 
address the problem of electronic structure and magnetic interactions in 
TMO on a more superior level is certainly a very important direction for 
the future. 

However, even in the present form, the electronic structure calcula- 
tions may have a significant impact on our understanding of the physical 
properties of TMO. In the next section we would like to illustrate how this, 
sometimes rather limited information extracted from the conventional elec- 
tronic structure calculations can be used for the analysis of many important 
questions accumulated in the field magnetism of TMO. 

5 Applications for Transition-Metal Oxides 

5.1 Double Exchange Interactions in CMR Manganites 

One of the most famous groups of oxide materials, which was in the fo- 
cus of enormous experimental and theoretical attention for many years, is 
perovskite manganese oxides (or simply - the manganites) with the chem- 
ical formulas Ri^xD^MnO^ (cubic type), Ri^xDi+xMnO^ (layered type), 
and i?2-2a;^i+22;Mn207 (doublc-laycrcd type; R and D being trivalent rare- 
earth and divalent alkaline-earth ions, respectively). The manganites are 
known for the colossal magnetoresistance effect, that is gigantic suppres- 
sion of the resistivity by the external magnetic field. Another important 
feature of these compounds is the rich magnetic phase diagram (Fig. HI) 
Typically, the CMR effect occurs at the boundary of either FM and PM, 
or FM and AFM phases. In this sense, the understanding of the magnetic 
phase diagram is directly related with the understanding of the CMR ef- 
fect. There is a large number of modern review articles devoted to these 
compounds [2Sl ESI EH E5|, which covers a lot of experimental informa- 
tion as well as the theoretical views on the problem. In this section we 
would like to illustrate, as we believe, the main idea behind the magnetic 
phase diagram and the CMR effect itself from the viewpoint of electronic 
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Figure 4: Magnetic phase diagram of Ndi_2,.Sra;Mn03 (from Ref. |S21)- Different phases 
are denoted as follows: PM - paramagnetic, FM - ferromagnetic, and AFM - antiferro- 
magnetic, consisting of the FM layers (A), simple chains (C), and zigzag chains (CE). 

structure calculations and the general theory of interatomic magnetic in- 
teractions presented in Sec. El We will concentrate on the doping range 

According to the formal valence arguments, the 3d states of Mn in 
Ri-xDxMnO^ are filled by the (4— x) electrons, which are subjected to the 
strong Hund's rule coupling. Therefore, the intra-atomic exchange split- 
ting between the f- and j-spin 3d states can be regarded as the largest 
parameter in the problem. The 3d states are further split by the cubic 
crystal-field into the lower-lying t2g and higher-lying Cg levels, yielding the 
formal electronic configuration t2gi^e^g^^t2gj^e^gi. The band dispersion, caused 
by the interatomic hopping interactions involving t2g and Cg orbitals is typ- 
ically smaller than the intra-atomic exchange and crystal-field splittings. 
The hoppings are mediated by the oxygen 2p states, and in the undistorted 
cubic lattice are allowed only between orbitals of the same, either t2g or e^,, 
symmetry. The strengths of the Mn{eg)-0{2p) and Mn{t2g)-0{2p) interac- 
tions are controlled by the Slater-Koster parameters, correspondingly pda 
and pdTT. Since |j9(ia'|>|j9(i7r|, the t2g band is typically narrower than the Cg 
one [36j. 

^^The magnetic phase diagram at small x is significantly affected by the Jahn- Teller effect. This is 
certainly an interesting and not completely understood problem (see, e.g., Ref. |26|V However, the physics 
is rather different from a more canonical double exchange regime realized for a;>0.3. 
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Thus, the t2g states will be half-filled and according to Sec. 13.11 con- 
tribute only to the AFM SE interactions. The Cg states are partially occu- 
pied and therefore will contribute to both FM DE and AFM SE interac- 
tions, whose ratio will depend on x. 

What is the relevance of this picture to the magnetic phase diagram 
of perovskite manganites? In order to get a very rough idea about this 
problem, let us simply count the number of FM (zfm) and AFM (^afm) 
bonds formed by the Mn sites in each of the magnetic phases shown in 
Fig. m (Table EI). The FM interactions prevail at smaller x. However, fur- 
Table 2: Number of ferromagnetic (zfm) and antiferromagnetic (-Zafm) bonds for the main 
magnetic phases observed in Ndi-^^Sr^^MnOa. 



phase 


X in phase diagram 


zfm 


zafm 


FM 


0.3<x<0.5 


6 





CE 


x~0.5 


2 


4 


A 


0.5<x<0.6 


4 


2 


C 


a;>0.6 


2 


4 


G 
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ther increase of x leads to the gradual change of the character of magnetic 
interactions towards the antiferromagnetism, which is refiected in the in- 
crease of the number of AFM bonds. Therefore, it seems reasonable to link 
this phase diagram with some kind of competition between FM and AFM 
interactions. 

In order to further proceed with this picture, we need to address the 
following questions. 

1. Which mechanism stabilizes the anisotropic AFM A and C structures 
in the metallic regime?^^ Note, that according to the canonical DE 
model PI , the metallic antiferromagnetism is unstable with respect to 
a spin-canting. 

^^According to the experimental data, the A phase is metaUic while the C phase is an insulator due to 
the cation disorder in the quasi-one-dimensional C-type AFM structure. Our point of view is that this 
disorder is not directly related with the magnetic stability of the C phase, and as the first approximation 
this phase can be treated in the conventional DE model which corresponds to the metallic behavior. 
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2. What is so special about the CE phase? Why does it exist only in 
the narrow window of x close to x=0.5 and perturbs the monotonous 
change of ^fm and ^afm shown in Table EI? Why is it insulator? 

3. Quantitative description of the doping-dependance of the magnetic 
phase diagram. 

5.1.1 Basic Details of the Electronic Structure 

The key factor factor, which allows to resolve the questions summarized in 
parts 1 and 2 is the strong dependance of the electronic structure of the Cg 
states on the magnetic structure EHl P9j . The effect can be estimated 
in the LSDA by considering the distribution of the Cg states in the hypo- 
thetical virtual-crystal alloy Lai/2Bai/2Mn03 (Fig. (5]). In this calculations, 
the crystal structure of Lai/2Bai/2Mn03 was rigidly fixed to be the cubic 
one. Hence, Fig. [51 shows the pure response of the electronic structure 
of the Cg states to the change of the magnetic structure. Even in LSDA, 
which typically underestimates the intra-atomic splitting between the f- 
and |-spin states jlOI , the effect is very strong and the electronic structure 
changes dramatically upon switching between FM, A-, and C-type AFM 
phases. Furthermore, the anisotropic AFM ordering has different effect on 
different Cg orbitals. For example, in the A-type AFM structure, the main 
lobes of the x'^-y'^ and orbitals are aligned correspondingly along FM 

and AFM bonds. Since the hopping interactions are penalized in the case 
of the AFM coupling, the A-type AFM arrangement will shrink mainly 
the 3z^-r^ band. The effect is reversed in the C-type AFM state, which is 
associated with the narrowing of the x'^-y'^ band. Thus, the LSDA calcula- 
tions clearly demonstrate how the anisotropy of the A- and C-type AFM 
structures interplays with the anisotropy of the spacial distribution of the 
x'^-y'^ and Sz'^-r'^ orbitals. 

The AFM CE phase presents the most striking example of such in- 
terplay between spin and orbital degrees of freedom: the zigzag AFM ar- 
rangement appears to be sufficient to explain the insulating behavior of 
this phase (even in LSDA). Therefore, the CE phase can be regarded as a 
band insulator, which is related with the very peculiar form of the AFM 
spin ordering (Fig. Ej) [SI- 
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Figure 5: Distribution of the Mn{eg) states in the ferromagnetic, A-, and C-type antifer- 
romagnetic phases of Lai/2Bai/2Mn03 in LSDA. The Fermi level is at zero. 

The magnetic-state dependence of the electronic structure of the Cg 
states is the main microscopic mechanism, which stands behind the rich 
variety of the magnetic structures observed in perovskite manganites. In 
some sense, these systems show certain tendency to the self-organization, 
the main idea of which can be accumulated in the following formula: the 
anisotropic AFM ordering breaks the cubic symmetry of the crystal and 
leads to the strong anisotropy of the electronic structure. This anisotropy 
will be reflected in the anisotropy of interatomic magnetic interactions (EI), 
which depend on the electronic structure through the matrix elements of 
the Green function. In many cases, the response of interatomic magnetic 

^■^Experimentally, this anisotropy is frequently observed as an orbital ordering 42 . 
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Figure 6: Right panel. Two-dimensional zigzag antiferromagnetic ordering realized in 
Lai/2Sr3/2Mn04, resulting in two different sublattices of the Mn sites (A and B, typically 
referred to the "charge ordering") and the orbital ordering (the alternation of partially 
populated 3x^-r^ and Sy'^-r'^ orbitals at the A sites). Left panel. Corresponding distribu- 
tion of the Mn(eg) states at the sites A and B in LSDA. The Fermi level is at zero. 

interactions to the cliange of the magnetic structure appears to be sufficient 
to explain the local stability of this magnetic state, and the total energy 
of such systems will have many local minima corresponding to different 
magnetic phases. 

5.1.2 Minimal Model for CMR Manganites 

The semi-quantitative model for CMR manganites can be constructed by 
using the following arguments. 

• The model should be based on the reahstic electronic structure for 
the eg states and take into account strong dependence of this elec- 
tronic structure on the magnetic structure. The proper model for 
the electronic structure can be obtained by doing the tight-binding 
parametrization of LDA bands 1201 SSI- However, for many applica- 
tions one can use even cruder approximation and consider only the 
NN hoppings between Mn sites, parameterized according to Slater 
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and Koster 144 



to to f . 27Ta 27Ta \ 

Ur+n^ " 2" ~ 2" \ ^ I ' 

where tT-T-'=||i:^^/ 1| is the 2x2 matrix in the basis of two Cg orbitals x^- 
y"^ (L,L'=1) and Sz^-r^ (L,L'=2), and is the primitive translation 
in the cubic lattice parallel to the x (a=l), y (a=2), or z {a=3) 
axis. The anisotropy of the x'^-y'^ and 3z^-r^ orbitals is reflected in 
the anisotropy of transfer interactions. The parameter to(x{pda)'^ is 
chosen to reproduce the Cg bandwidth in LDA (to— 0.7 eV [i26j, and 
thereafter used as the energy unit). 

The effect of the magnetic structure on the electronic structure is 
typically included after transformation to the local coordinate frame, 
specified by the directions {gt-} of the spin magnetic moments 

Gt = (cos (j)-r sin sin (f)^ sin 6^-, cos 6^-) , 

and taking the limit of infinite intra-atomic exchange splitting. This 
yields the well-known DE Hamiltonian [45j : 

H^^t = —^TT'tTT', (12) 

where 

F^^i = COS — COS h sm — sm — e ' 

^ 2 2 2 2 

describes modulations of transfer interactions caused by deviations 

from the FM spin alignment. 

• 'H^-t' proportional to to, and similar to the DE term considered in 
Sec. 13. IL describes the FM interactions in the systems, which are now 
generalized to the case of arbitrary spin arrangement. This term is 
combined with the energy of phenomenological AFM SE interactions: 



^'^The DE Hamiltonian is essentially non-local with respect to the site indices. In this sense one can say 
that the DE physics is mainly the physics of bonds, which may have a direct implication to the behavior 



of paramagnetic phase of CMR manganites |46|. 
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where J^<0 is restricted by the nearest neighbors and contains con- 
tributions of both t2g and Cg states. Using an analogy with the SE 
model considered in Sec. 13. 21 it is assumed that the exchange constant 
does not depend on the magnetic state, and all such dependencies 
come exclusively from the DE term. 

The combination of H^^f and H^^f constitutes the minimal model, 
which explains the origin of the main magnetic structures observed in 
doped manganites.^^ An example of theoretical phase diagram is shown in 
Fig. [71 |37l EBl EDI • The position and the order of the main magnetic phases 
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Figure 7: The areas of local stability of the main coUinear phases with respect to the spin 
canting in the degenerate double exchange model (from Ref. 08:). 

with respect to the hole-doping has clear similarity with the experimental 
phase diagram (Fig. HI). 

The origin of the CE-type AFM phase requires special comment, be- 
cause canonically it was attributed to the charge ordering of the Mn^+ and 
Mn^"*" ions [Ij, and this point of view was dominating for many years. The 
modern view on the problem is that CE is the regular magnetic phase, 
whose properties can be understood in the degenerate DE model similar 

^''The on-site Coulomb interaction |17] and the lattice distortion 021 have been also considered as 
the model ingredients. However, as far as the low-temperature behavior is concerned, the inclusion of 
these terms (of course, under the appropriate choice of the parameters 0^1) may change the situation 
quantitatively, but not qualitatively. 

i^Note that around x=0.5, | J'^| is of the order of 0.1 EOl- 
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to other magnetic states [2Sl EHl EHj • The insulating behavior of this phase 
is related with the topology of one-dimensional zigzag FM chains. This 
also explains why the CE phase exists in the very narrow region close to 
x=0.5, i.e. when the Fermi energy falls into the band gap. More generally, 
any periodic one-dimensional zigzag object composed of the Cg orbitals will 
behave as a band insulator at certain commensurate values of the electron 
density [2BIEI]. 

5.1.3 Implication to the CMR Effect 

The magnetic origin of the CE phase has many implications to the CMR 
effect and readily explains why the insulating state of some manganites can 
be easily destroyed by the external magnetic field. However, it would not 
be right to say that CMR is exclusively related with the existence of the 
CE phase. It is a more general phenomenon inherent to the DE physics. 
Below we will consider two such possibility. 

1. According to the form of the DE Hamiltonian (IT^. the transfer inter- 
actions are forbidden between antiferromagnetically coupled sites. There- 
fore, from the viewpoint of transport properties, the A-type AFM phase 
behaves as an insulator in the z-directions and as a metal in the xy-plane. 
The A state can be continuously transformed to the spin-canted state |2B] 
by applying the external magnetic field. The spin-canting unblocks the 
hopping interactions and gradually decreases the resistivity along z. This 
is the basic idea of the so-called "spin- valve" effect, observed in the A-type 
AFM Ndo.45Sro.55Mn03 |S2|. 

The unique aspect of the CE phase is that it has a band gap, and there- 
fore is insulating in all the directions, including the direction of propagation 
of the FM zigzag chains. This predetermines rather distinguished behavior 
of the CE phase in the magnetic field. As it was pointed out before, the 
stability of the CE phase is directly related with the existence of the band 
gap. Similar to the A-type AFM phase, the external magnetic field leads 
to the canting of spins also in the CE phase. The basic difference is that 

^^The Cg electron passing through the 3z^-r^ and x^-y^ orbitals of the corner (B) sites (Fig. O cor- 
respondingly preserves and changes its phase. In the one-dimensional case, it will open the band gap 
correspondingly at the boundary and in the center of the Brillouin zone, that explains the insulating 
behavior. 
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starting with certain angle such canting wiU close the band gap. Then, 
the CE phase will become unstable and the system will abruptly transform 
to the metallic FM phase, which is favored by the magnetic field. The phe- 
nomenon is well know as the " melting of the charge-ordered state" , which 
accompanied by the abrupt drop of the resistivity [34J . 

2. Typical theories of the second type are based on the following ob- 
servations fH^:^^ 

• The total energy of the DE model may have several local minima 
corresponding to different FM and AFM phases, which may exist in 
the same range of the hole-doping. This also implies the first-order 
character of the phase boundaries jSD|- The concrete example is the 
region close to x=0.5 shown in Fig. [71 

• Some of these phases are insulating due to the special (zigzag) geom- 
etry of the AFM pattern. 

The situation is schematically illustrated in Fig. |Hl At T=0 and without 
field, the system is trapped in the insulating CE-type AFM ground state. 
The next minimum, corresponding to the FM metallic state should be 
located within the energy range accessible by the external magnetic field. 
At the finite temperature and neglecting interactions between different 
phases, the system will be represented by the mixture of these two states, 
due to the configuration mixture entropy. The relative amount of these 
states will depend on the temperature T and the magnetic field H. The 
latter controls the relative position of the CE and FM minima. Without 
field, the CE phase will dominate, while the FM phase will form metallic 
islands in the insulating sea. The application of the magnetic field will 
increase the amount of the of the FM phase. At some point, which is called 
the percolation threshold, the FM islands become connected by forming the 

-"^^In the two-dimensional DE model, the band gap is closed when the angle between spin magnetic 
moments in adjacent zigzag chains becomes smaller than 110°. 

^^The percolative scenario of the CMR effect in manganites has been originally proposed by Nagaev )54| . 
The newer theories provide a quantitative description of this effect and clarify the origin of the main 
magnetic states, particularly - the existence of the insulating CE-type AFM phase in the degenerate DE 
model. 

^°The appearance of the two-phase state can be considerable facilitated in the presence of impurities 
and the cation disorder [531 154) . As we will argue in Sec. 15.21 this effect plays a very important role in 
the double perovskites. 



26 




CE 




CE 


FM 


1 

^ /fm 


CE 






CE 


J FM \ 


f FM ^ 

1 


/ C^E 



Figure 8: Schematic view on the percolative scenario of the CMR effect. Top: Local 
minima of the double exchange energy corresponding to the insulating CE and metallic 
FM phases. The relative position of these minima is controlled by the magnetic field 
H. Bottom: mixed-phase state realized at the finite temperature. The magnetic field 
increases the amount of the FM phase, and above the percolation threshold forms the 
conducting FM path. 

conducting FM paths throughout the spacemen. This wiU correspond to 
the sharp drop of the resistivity. 

Since many interesting (and, presumably, the most practical) phenom- 
ena of perovskite manganites are related with the existence of zigzag AFM 
structures, it is very important to understand the behavior of these struc- 
tures, and especially the way how they will evolve with the temperature, in 
the magnetic field, or upon the change of the crystal structure. We would 
like to close this section by listing some interesting and not completely 
resolved problems in this direction. 

• The existence of two transition temperatures, Tn and Tco, which are 
typically attributed to the onset of Neel AFM and charge ordering, 
respectively. In reality, both transitions are probably magnetic and 
Tn corresponds to the order-disorder transition which takes place only 
between FM zigzag chains and largely preserves the coherent motion 
of spins within the chains, while Tco corresponds to the transition to 
the totally disordered FM state 



27 



• Appearance of incommensurate charge and orbital superstructures 
just below Tco IBS]- 

• Effects of cation disorder and appearance of new zigzag superlattices 
upon artificial ordering of R and D elements in some three-dimensional 
manganites [H7j . 

5.1.4 Role of Rare-earth, Alkaline-earth, and Oxygen and States 

In this section we will briefly discuss the influence of the i?, and oxy- 
gen states on the interatomic magnetic interactions in Ri_xDxMnO^ in the 
context of revisions in the form of these interactions, which were recently 
considered by Bruno [201 • We will concentrate on the behavior of the cu- 
bic FM phase of the virtual-crystal alloy LaojBao.sMnOa, and treat the 
problem of interatomic magnetic interactions using the frozen spin wave 
approximation supplemented with the LSDA.^^ In the FM phase, the mag- 
netic Mn sites polarize neighboring La/Ba and oxygen sites. The typical 
distribution of the magnetic moments amongst different sites is /iMn=3.48, 
MLa/Ba=0-08, and /io=0.05 fiB- Our main concern will be the role played 
by /tiLa/Ba Siud fiQ in the spin dynamics of Lai-^Ba^^MnOs. 

We assume that the cone-angle of the spin wave does not depend on 
the atomic site: 9t-=6. The phase of the spin wave is modulated as q • r 
with T running over all Mn, O, and La/Ba sites. Then, for each q we 
calculate the second derivative of the total energy with respect to 6 using 
the magnetic force theorem. The derivative can be also used to estimate 
the magnon spectrum in the FM state: 

_ 2 £E{q,e) 

{fi being the FM moment). We consider two possibility: 

• Rigid constraint, obtained by rotating the KS effective field and en- 
forcing the directions of the spin magnetic moments at each site to 
follow exactly the form of the spin wave. The latter is achieved by ap- 
plying the perpendicular constraining fields {h:;^} (Fig. IHI), which can 

^^The cubic lattice parameter is ao=3.876A. 
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Figure 9: Typical situation realized in frozen spin wave calculations for q=(0, 0, vr). {/x} is 
the constrained distribution of the magnetic moments, {/i'} is the distribution of magnetic 
moments obtained after solution of Kohn-Sham equations with the rotated effective fields. 
Generally, the spins {a*'} are canted towards the ground state configuration by the angle 
T), and have the components {a*"*"}, which are perpendicular to {/x}. The constraining 
fields {h^} are included in order to compensate {m^}. 

be estimated using the response-type arguments [2D]- The fields {h:;^} 
affect the KS orbitals {V^i}, and through the change of these orbitals, 
correct the form of the spin-magnetization density (I2j). However, they 
do not contribute explicitly to the total energy change (EI). 

• Relaxed constraint, obtained after rotation of the KS effective fields. 

Results of these calculations are shown in Fig. [TUl The magnon dispersion 
has a characteristic form, which is manifested in the pronounced softening 
at the Brillouin zone boundaries [58j.^^ The origin of this softening is 
related with the behavior of interatomic magnetic interactions, and can be 

^^Experimentally, the softening of the spin- wave dispersion is frequently accompanied by a broadening, 
though these two effects are not necessarily related with each other and the materials which show very 
similar softening of the magnon spectrum can reveal very different broadening |58| . 
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Figure 10: Results of frozen spin wave calculations for Lao.rBao.aMnOs: second deriva- 
tives of the total energy (which specify the form of the magnon dispersion) obtained by 
rotating the Kohn-Sham effective fields (relaxed constraint) and by applying additional 
perpendicular constraining fields {h;^} (rigid constraint) at all sites and the Mn sites only; 
the parasitic canting of the magnetic moments in the relaxed constraint calculations; and 
the perpendicular magnetic field applied in order to compensate this canting. Symbols 
show the calculated points, and the solid lines are results of an interpolation. 

understood on the basis of realistic electronic structure calculations [25j. 
It has a direct implication to the stability of the FM state and the form of 
the magnetic phase diagram versus the hole-doping jlHI • 

In this section, our main concern will be a little bit different. As it is 
seen in Fig. [TUl the application of the constraining fields {h:;^} can signif- 
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icantly affect tfie magnon dispersion. For example, around (tt, 0, 0), cjq's 
calculated using two different schemes can differ by factor two. The most 
interesting question is where this difference comes from. The simple rota- 
tion of the KS effective fields by the angle 9 causes the parasitic canting 
of spins off the constrained directions by the angles {r/r}, which are pro- 
portional to 9 if the latter is small (see Fig. IHI). However, this canting can 
be very different at different magnetic sites. Indeed, rjr is of the order of 
9 at the La/Ba and oxygen sites, and almost negligible at the Mn sites. 
Therefore, the fields {h:;^}, which are introduced in order to correct this 
canting will be the largest at the La/Ba and oxygen sites, and the smallest 
at the Mn sites. An attempt to apply h:;^ only at the Mn sites has a negli- 
gible effect on the magnon dispersion. These results rise several important 
questions, which still need to be understood. 

• There is a large intrinsic uncertainty in the behavior of spin magneti- 
zation at the La/Ba and oxygen sites, and especially in the form, in which 
these sites contribute to the spin dynamics of Lai_a;Ba2;Mn03 and similar 
compounds. The problem may be directly related with the definition of 
the adiabaticity concept, underlying the frozen spin wave approach. Ac- 
cording to this concept, the magnetic systems have two characteristic time 
scales. The spin directions are typically regarded as the slow variables, 
while the internal electronic degrees of freedom are the fast ones, which 
self-consistently follow the instantaneous orientational distribution of the 
spin-magnetization density. From this point of view, the La/Ba and oxygen 
sites have clearly a dual character. On the one hand, neither of them can be 
regarded as the source of magnetism in these systems. On the other hand, 
both of them carry a significant weight of the spin-magnetization density, 
mainly due to the hybridization with the magnetic Mn sites. Then, how 
the La/Ba and oxygen states shall be treated? In principle, one may have 
several possibihties. 

1. The La/Ba and oxygen sites are magnetic and shell be treated on 
the same level as the Mn sites, by applying independent constraint 
conditions for the slow- varying directions of the La/Ba and oxygen 
spins. This would correspond to the appearance of additional magnon 
branches in the frozen spin wave calculations. Clearly, this scenario 
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can be sorted out from the beginning, because both La/Ba and oxy- 
gen moments are induced by the hybridization with the Mn sites and 
cannot be totally independent from the latter. 

2. The contribution of La/Ba and oxygen states is purely electronic, 
which implies that these degrees of freedom are "fast" and self- 
consistently follow the distribution of the spin magnetic moments at 
the Mn sites. This would roughly correspond to the relaxed constraint 
conditions in the spin- wave calculations shown in Fig. [TDl 

3. An intermediate situation when the electronic degrees of freedom at 
the La/Ba and oxygen sites are not "sufficiently fast" so that the 
magnon dispersion shall be averaged over several intermediate config- 
urations of the La/Ba and oxygen spins. This is the most interesting 
scenario, which is however only a speculation. If true, the magnon dis- 
persion in FM manganites should have an intrinsic bandwidth, which 
should be of the same order of magnitude as the difference between 
"rigid" and "relaxed" calculations shown in Fig. [mi In fact, the strong 
magnon broadening has been observed experimentally in several FM 
manganites [SB|, though it is of course premature to make a direct 
connection between this broadening and the behavior of La/Ba and 
oxygen states considered in this section. 

• The non-collinear distribution of magnetic moments at the La/Ba and 
oxygen sites can interfere with other (extrinsic) factors, such as randomness 
effects caused by the cation disorder. Then, the chemical disorder may 
lead to the orientational spin disorder at the La/Ba and oxygen sites, 
which may also contribute to the magnon broadening associated with the 
randomness [HSj. 

5.2 Double Perovskites 

Ordered double perovskites, like Sr2FeMo06 and Sr2FeRe06, is another 
interesting class of compounds [60j. They exhibit fairly large magneto- 
resistance effect, which coexists with relatively high magnetic transition 
temperature (415 and 401 K for Sr2FeMo06 and Sr2FeRe06, respectively). 
The latter factor presents a great advantage over CMR manganites, and 
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opens a way to technological applications of these compounds in devices 
operating at the room temperature. 

At the first sight, the double perovskites crystallize in a simple crystal 
structure. It can be regarded as the cubic SrFeOs, in which every second 
Fe is replaced by Mo or Re (Fig. fTTl). However, there are two fundamental 
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Figure 11: Positions of Fe, Mo, and oxygen sites in ordered double perovskite Sr2FeMo06. 
Arrows show the main magnetic interactions responsible for properties of this compound. 

problems: 

1. Anti-site defects. Typically, certain percent of atoms from the nomi- 
nally Fe sublattice interchanges with the same amount of atoms from 
the Mo sublattice. The number of such defects in the best available 
single-crystalline samples of Sr2FeMo06 is about 8% jH] . 

2. Fine details of the crystal structure. Typically, the information about 
directions and magnitudes of the oxygen displacement is either un- 
known or rather controversial. 

The extraordinary properties of Sr2FeMo06 and Sr2FeRe06 are usually 
attributed to the half-metallic electronic structure expected in the FM 
state. In this sense, the double perovskites are typically considered as an 
example of relatively simple and well established physics. 

In this section, we will discuss the behavior of interatomic magnetic in- 
teractions in Sr2FeMo06 and argue that this point of view must be largely 
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revised. Particularly, we will show that the half-metallic electronic struc- 
ture is incompatible with stability of the FM ground state. Therefore, 
there should be some additional mechanism which stabilizes the FM state 
and presumably destroys the half-metallic character of the electronic struc- 
ture. We propose that the unique properties of the double perovskites are 
directly related with the above-mentioned defects and uncertainties of the 
crystal structure |^.^^ 

A rough idea about the problem can be obtained from rather simple 
considerations of the electronic structure of Sr2FeMo06 (Fig. [T2I), which 
suggests that similar to the CMR manganites there are two competing 
interactions which define the form of the magnetic ground state in this 
compound. 

1. It is true that the electronic structure of the hypothetical FM phase of 
Sr2FeMo06 is half-metallic: the Fermi level crosses the hybrid t2g band 
in the j-spin channel and falls in the gap between Fe(eg) and Mo{t2g) 
bands in the j-spin channel. The partly-filled |-spin t2g band will gen- 
erally be the source of FM interactions, due to the double exchange 
mechanism considered in Sec. 13. IL which will be somewhat modified by 
the fact that the transfer interactions between spin-polarized Fe{t2g) 
orbitals are indirect and operate via nonmagnetic Mo{t2g) states [63j. 
This interaction will also polarizes the Mo(t2g) states antiferromag- 
netically with respect to the Fe sublattice. 

2. The half-filled Fe(e^) states should generally contribute to AFM in- 
teractions. Typically, this contribution is expected to be small, be- 
cause in the double-perovskite structure, the neighboring Fe sites are 
separated by the long 0-Mo-O paths. A very rough idea about the 
strength of these interactions can be obtained from the comparison of 
densities of the Fe{eg) states in the FM and (type-II [6j) AFM phases: 
if these interaction were weak, the bandwidth would not depend on 
the magnetic state. However, the direct calculations show the oppo- 
site trend and the Fe(e^) band considerably shrinks in the case of the 
AFM alignment. 

^•^The same arguments can be applied to Sr2FeRe06- 
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Figure 12: Local densities of states corresponding to the ferromagnetic (left) antiferro- 
magnetic (right) alignment of Fe spins in ordered Sr2FeMo06, in LSDA. The Fermi level 
is at zero. 

A simple estimate of the relative strength of the FM and AFM contri- 
butions, irrespective on their decomposition onto the interatomic magnetic 



35 



interactions, can be obtained from the one-electron (band) energies 



{J\f being the total density of states) calculated separately for the Cg and t2g 
bands in the FM and AFM phases. These energies are listed in Table El^^ 
One can clearly see that in the AFM state, the energy loss associated with 

Table 3: One-electron energies associated with the hybrid t2g band located near the Fermi 
level and the occupied Fe(eg) band in the FM and AFM states of Sr2FeMo06. 



the hybrid t2g band is totally compensated by the energy gain associated 
with the Fe(eg) band. This example shows that the problem is indeed 
very subtle and the half-metallic electronic structure obtained in the FM 
state does not necessary guarantee that this state will be locally stable 
and realized as the ground state. It does not explain the high value of the 
magnetic transition temperature either. 

The conclusion is supported by direct calculations of interatomic mag- 
netic interactions using both frozen spin wave and Green's function meth- 
ods. 

1. The FM state appears to be unstable with respect to the 
non-coUinear spin-spiral alignment realized at certain q's, for which 

E{q, O)\0^Q<O (Fig. [131). The second derivative of the total energy was 
calculated using several different approximations, similar to the analysis 
of the spin-wave dispersion of LaojSro.sMnOa undertaken in Sec. 15.1.41 by 
imposing the rigid constraint on the directions of spin magnetic moments 
at all sites of Sr2FeMo06, and only at the Fe and Mo ones. Similar to 
the LaojSro.sMnOa case, the main effect of the constraining fields {h.^} is 
associated with the Sr and oxygen sites, whereas the contributions of the 
Fe and Mo sites are negligible. The "rigid constraint" can be regarded 

^''Note that in the case of SE interactions, in addition to the change of the Fe(eg) band energy, there 
wiU be an additional contribution associated with the shift of the 0{2p) states fSec. 13.21 and Ref. 
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Figure 13: Results of frozen spin wave calculations for ordered Sr2FeMo06 (see Fig. HUl 
for details). 



as the most consistent check for the local stability of the FM state, when 
the application of the magnetic force theorem is strictly justified and the 
second derivative of the total energy can be rigorously expressed through 
the change of the KS single-particle energies. Fig. [T31 shows that it may 
correct to a certain extent the values of the second derivatives at the Bril- 
louin zone boundary. However, it does not change the principal conclusion 
and the FM state remains unstable in the large portion of the Brillouin 
zone 102!.^^ 

2. The second derivatives of the total energy can be mapped onto the 
Heisenberg model ^ £'(q, ^)|g,^Q=Jo— Jq. The mapping gives us parame- 
ters of magnetic interactions in the reciprocal space {Jq}, which can be 
further transformed to the real space. The real-space parameters, listed in 
Table HI provide a complementary piece of information which allows to ra- 
tionalize results of frozen spin wave calculations and present them in terms 
of competition between FM NN and AFM next-NN interactions in the Fe 
sublattice (correspondingly Ji and J2 in Table HI). The AFM interaction 
J2 clearly dominates and readily explains why the FM state becomes un- 
stable. An additional piece of information can be obtained from Green's 

^^Note that in order to prove that the magnetic state is unstable, it is sufficient to find at least one 
configuration of the spin-magnetization density for which ^i?(q, 9)<0. 

^^Note also that Ji and J2 have different coordination numbers (correspondingly 12 and 6). 
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Table 4: Parameters of magnetic interactions in the real space obtained in the (relaxed) 
frozen spin wave (SW) and Green's function (GF) calculations near the FM and AFM 
states of Sr2FeMo06 (in meV). Two values of Ji in the case of the AFM alignment corre- 
spond to nearest-neighbor interactions between Fe sites having the same (the first number) 
and opposite (the second number) spins in the type-II antiferromagnetic structure. 



method 


Ji 


J2 Jpe-Mo 


SW (FM) 


9.3 


-26.9 


GF (FM) 


6.2 


-21.5 -1.2 


SW (AFM) 


16.4, 18.6 


-15.2 



function calculations, which allow to estimate the direct Fe-Mo interaction 
Jpe-Mo- However, this interaction is not particularly strong and does not 
alter the main conclusion about stability of the FM state. 

3. If the FM state is unstable, what is the true magnetic ground state 
of the ordered Sr2FeMo06? The question is not simple, because the AFM 
phase appears to be also unstable, as it is expected for small concentrations 
of {t2g) carriers moving in the AFM background of localized (e^) spins and 
interacting with the latter via the strong intra-atomic Hund's coupling [7J. 
This fact is also related with the magnetic-state dependence of interatomic 
magnetic interactions, which leads to the reversed inequality Ji>|J2| for 
parameters calculated in the AFM state (Table HI). Therefore, the true 
ground state of the chemically ordered Sr2FeMo06 should he in-between the 
FM and AFM states. The character of interatomic magnetic interactions 
suggests that probably it is a spin-spiral with q||[lll]. 

The FM ordering can be stabilized by the anti-site defects or the local 
lattice distortions |iQ2j- However, both mechanisms will destroy the half- 
metallic character of the electronic structure of Sr2FeMo06. This is an 
irony of the situation. From the electronic structure point of view, the 
most efficient way to stabilize the ferromagnetism is to create holes in the 
|-spin Fe(eg) band, by shifting the Fermi level to the lower energy part of 
the spectrum, and thereby activate a very effective channel for the FM DE 
interactions associated with the eg states, similar to the CMR manganites. 

^^Small difference of the parameters obtained in the frozen spin wave and Green's function approaches 
is related with the fact that Jfo-Mo is already included in the definition of Ji and J2 in the frozen spin 
wave calculations. 
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Let us give a rough estimate for this effect. The width of the Fe{eg) 
band (W) in Sr2FeMo06 is about 2 eV (Fig. [T2I), which is approximately 
two times smaller than the typical Cg bandwidth in FM manganites. There- 
fore, the DE and SE interactions in Sr2FeMo06, being proportional to W 
and W"^, will be correspondingly two and four times smaller in comparison 
with the same interactions in manganites. The latter can be estimated 
as ~ 80 and ~ —70 meV, respectively PB|. Therefore, the partial 
depopulation of the |-spin Fe(e^) band in Sr2FeMo06 leads the following 
estimate for next-nearest-neighbor coupling: J2 — | J'^ + ^ J*^ ~ 20 meV, 
which easily explains the appearance of ferromagnetism in this compound. 

Thus, the half-metallic electronic structure realized in the FM state 
of the ordered double perovskites is a spurious effect since it is largely 
incompatible with stability of this FM state. It also implies that the 100% 
ordered double perovskites cannot be ferromagnetic. The latter observation 
is qualitatively consistent with the behavior of Sr2FeW06, where almost 
perfect ordering of Fe and W coexists with an AFM ground state [.64J.^^ 

The quantitative theory of ferromagnetism in Sr2FeMo06 and 
Sr2FeRe06 is still missing. The exceptionally high value of the Curie tem- 
perature is one of the main puzzles, which is probably related with the 
extrinsic inhomogeneities (anti-site defects, local lattice distortions, grain 
boundaries) existing in the samples. The magneto-resistive behavior is 
based on the percolative scenario, and is not necessary related with the 
half-metallic character of the electronic structure. The main difference 
from CMR manganites is that neither FM nor AFM states are locally sta- 
ble and cannot be the total energy minima in the case of double perovskites. 
Therefore, the extrinsic inhomogeneity seems to be indispensable in order 
to stabilize the FM islands and form the two-pase state similar to the one 
depicted in Fig. |Hl^^ Then, the negative magneto-resistance may occur if 
the size of these islands can be easily controlled by the external magnetic 
field. 

Finally, our analysis was based on the electronic structure obtained in 

^®The assignment is based on the analysis of magnetization data, which cannot exclude the non-coUinear 
spin-spiral alignment. 

^^This was also supported by experimental studies of Sr2FeWi_2;Moa;06 alloy, which shows clear ten- 
dency to the phase segregation )64| . The largest magneto-resistance was observed around a;=0.15. 
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LSDA, which may be imperfect. The role of Coulomb correlations on the 
top of LSDA is rather unclear. It is also unclear whether the Coulomb 
correlations alone can stabilize the FM state. On the one hand, the on-site 
U enters the denominator for the SE interactions and therefore suppress the 
AFM contribution associated with the Fe(eg) states. On the other hand, 
it may change the relative position of the Fe and Mo states and thereby 
suppress the FM DE interactions associated with the hybrid t2g band [021 • 
We expect that the LSDA band structure depicts the main problems of 
the double perovskites, though the concrete details may depend on the 
Coulomb correlations on the top of this picture. 

5.3 ^aMoaOy {A= Y, Gd, and Nd) Pyrochlores 

Pyrochlores with the chemical formula 742M02O7 (^4 being the divalent el- 
ement) is one of the rare examples of magnetic oxides on the basis of 4(i 
elements. The magnetic phase diagram of 742M02O7 is controlled by the 
averaged ionic radius {ta) of the 74-sites. This dependence is very puzzling 
(Fig. ED [661 EH Ea t69j. Large (r^) {>Rc^lMl) stabilizes the ferro- 
magnetic (FM) ground state. Typical examples of the FM pyrochlores are 
Nd2Mo207 and Gd2Mo207. The Curie temperature (Tc) is of the order of 
80 K and slowly increases with (r^). Smaller [r a) {<Rc) give rise to the 
spin-glass (SG) behavior. The characteristic transition temperature to the 
SG state is of the order of 20 K. The canonical example of the SG com- 
pounds is Y2M02O7 [TDj. Intuitively, it is clear that such a behavior should 
be related with the change of the NN exchange coupling, which becomes 
antiferromagnetic in the SG region. The pyrochlore lattice supplemented 
with the AFM exchange interactions presents a typical example of geomet- 
rically frustrated systems with an infinitely degenerate magnetic ground 
state [7T], which probably predetermines the SG behavior. The FM-SG 
transition in Mo pyrochlores is closely related with the metal-insulator 
transition. All compounds which reveal the SG behavior are small-gap 

Another famous example of magnetic 4d oxides is the ruthenates, crystallizing in the cubic, layered, 
and double-layered perovskite structures |65j . 

■^^Nevertheless, many details of this behavior remain to be understood. For example, according to the 
recent experimental data [72], the SG state appears due to the joint effect of geometrical frustrations 
and the disorder of local lattice distortions. The latter is responsible for the spacial modulations of 
interatomic magnetic interactions, which freeze the random spin configuration. 
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Figure 14: Experimental phase diagram of yl2Mo207: transition temperature against the 
averaged ionic radius. The dotted hne shows the boundary between normal ferromagnetic 
(right) and spin-glass (left) behavior. 

insulators. However, the opposite statement is generally incorrect and 
in different compounds the ferromagnetism is known to coexist with the 
metallic as well as the insulating behavior jHSj- In this section we will try 
to understand which parameter of the crystal structure controls the sign of 
the NN magnetic interactions, and which part of the electronic structure 
is responsible for the FM and AFM interactions in these compounds. 

5.3.1 Main Details of Crystal and Electronic Structure 

The pyrochlores 742M02O7 crystallize in a face-centered cubic structure 
with the space group Fm3d, in which A and Mo occupy correspondingly 
16d and 16c positions, and form interpenetrating sublattices of corner- 
sharing tetrahedra. There is only one internal parameter which may control 
the properties of 742M02O7. That is the coordinate u of the O 48/ sites. 

The single Mo tetrahedron is shown in Fig. [T51 Four Mo sites are 
located at Ti=(0,0,0), t2=(0, T3=(i,0,i), and t4=(|, |,0), in units 
of cubic lattice parameter a. Each Mo site has sixfold O 48/ coordination. 
The oxygen atoms specify the local coordinate frame around each Mo site. 
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Figure 15: Positions of the Mo and O 48/ sites in the pyrochlore lattice. 



Around 1, it is given by 

nf = (13) 

where a, (3= x, y, or z. u=jq corresponds to the perfect octahedral en- 
vironment, while u>jq gives rise to an additional trigonal contraction of 
the local coordinate frame. Similar matrices associated with the sites 2, 
3, and 4 can be obtained by the 180° rotations of TZi around x, and z, 
respectively. 

Structural parameters of 742M02O7 (taken from Refs. |^ 1^ IHHll^ 170] ) 
are listed in Table El Corresponding densities of states, obtained in LSDA 
are shown in Figs. [121 and [T71 

In the local coordinate frame, the Mo(4(i) orbitals are split into the 
triply-degenerate t2g and doubly-degenerate Cg states. Twelve t2g bands 
are located near the Fermi level and well separated from the rest of the 
spectrum. Interestingly enough that all three compounds 742M02O7 (^4= 
Y, Gd, and Nd) are ferromagnetic, even in LSDA, that is rather unusual for 
the M oxides, perhaps except the well know example of SrRuOs [73]. The 
trigonal distortion and the different hybridization with the 0{2p) states 
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Table 5: Structural parameters of A2M02O7 {A= Y, Nd, and Gd): the cubic lattice 
parameter a (in A), positions of O 48/ sites [u, |, |] (units of a), the distances Mo- Mo (in 
A), and the angles Mo-O-Mo (in degrees). 



compound 


a 


u 


C^Mo-Mo 


ZMo-O-Mo 


Y2M02O7 


10.21 


0.33821 


3.6098 


127.0 


GdsMosOy 


10.3356 


0.33158 


3.6542 


130.4 


NdsMoaOr 


10.4836 


0.32977 


3.7065 


131.4 



will further split the Mo{t2g) states into the one-dimensional aig and two- 
dimensional e'g representations.^^ The crystal structure affects the Mo{t2g) 
bandwidth via two mechanisms (see Table EI). 

1. The Mo-O-Mo angle, which decreases in the direction Nd^Gd^Y. 
Therefore, the interactions between Mo{t2g) orbitals which are medi- 
ated by the 0{2p) states will also decrease. 

2. The lattice parameter a and the Mo-Mo distance, which decrease in 
the direction Nd^Gd^Y by 2.6%. This will increase the direct in- 
teractions between extended Mo(4(i) orbitals. 

Generally, these two effects act in the opposite directions and partly com- 
pensate each other. For example, the width of the e'g band is practically 
the same for all three compounds (Fig. fTTl) . On the other hand, the aig 
orbitals, whose lobes are the most distant from all neighboring oxygen sites 
are mainly affected by the second mechanism, and the aig bandwidth will 
increase in the direction Nd^Gd^Y. As we will see below, this effect plays 
a crucial role in the stabilization of AFM interactions in Y2M02O7, which 
predetermines SG behavior. 

Thus, despite an apparent complexity of the crystal structure, the py- 
rochlores 742M02O7 present a rather simple example of the electronic struc- 
ture and in order to understand the nature of fascinating electronic and 
magnetic properties of these compounds, we need to concentrate on the 
behavior of twelve well isolated Mo{t2g) bands. 

■^^In the local coordinate frame, the aig and two e'g orbitals have the following form: 
\aig) = ^i\xy) + \yz) + \zx)), \e'gi)=^i\yz)-\zx)), and \e'g^) = ^{-2\xy) + \yz) + \zx)). 
^■^In this sense, the physics is similar to the spinel compounds [T^. Note, however, that the oxygen 
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Figure 16: Total and partial densities of states of Nd2Mo207 in the local-spin-density 
approximation. The Mo(t2g) states are located near the Fermi level (chosen as zero of 
energy). The Mo{eg) states emerge around 4 eV. 

In this section we consider results of model Hartree-Fock calcula- 
tions [7S], which combine fine details of the electronic structure for these 
bands, extracted from LSDA, and the on-site Coulomb interactions among 
the Mo(4(i) electrons, treated in the most general rotationally invariant 
form [I3|.^'^ 

coordination is very different in the spinel and pyrochlore structures. 

'^''The on-site Coulomb interactions are specified by three radial Slater integrals, F^, and F°. 
Equivalently, one can introduce parameters of the averaged Coulomb interaction, U —F'^, the intra-atomic 
exchange coupling, J=j^{F^+F'^), and the non-sphericity of Coulomb interactions between orbitals with 



44 



8 




-1,2-1.0-0,8-0.6-0.4-0,2 0,1) 0,2 0.4 0.6 0.8 1,0 
Energy (eV) 



Figure 17: Mo(t2g) states in the local coordinate frame, split into the one-dimensional aig 
and two-dimensional e'g=eg{t2g) representations by the local trigonal distortion. 

5.3.2 Effects of Crystal Structure and On-Site Coulomb Interactions on the 
Interatomic Exchange Coupling 

According to the LSDA calculations for the FM state (Fig.[T7I), the |-spin 
aig band is fully occupied and the Fermi level crosses the doubly-degenerate 
e'g band. Therefore, it is clear that at some point the Coulomb U will split 
the e'g band and form an insulating state with the spontaneously broken 
FrnSd symmetry. Such situation occurs between U= 2.0 and 2.5 eV for 
all considered compounds (Fig. [TBI). Typical densities of states in the 

the same spin, B=^{9F^—5F^). J~0.5 eV is taken from the constraint-LSDA calculations B 
can be estimated from J using the ratio F''^/F^~0.63, which holds for the Slater integrals in the atomic 
limit. This yields B~0.06 eV. The Coulomb U is treated as the parameter in order to consider different 
scenarios, covering both metallic and insulating behavior of A2M02O7. The constraint-LSDA calculations 
for Mo compounds yield J7«3.0 eV This value can be further reduced by allowing for the (proper) 
eg electrons to participate in the screening of on-site Coulomb interactions of the t2g electrons |7(i| . 



45 



> 



' 0.0 0.5 1,0 1.5 2.0 2.5 3,0 3.5 
i7(eV) 

Figure 18: The band gap as a function of Coulomb f/. For f/<2.0 eV there is an overlap 
between bands, corresponding to the negative value of the band gap, while f/>2.5 eV 
opens the real gap. 



insulating phase are shown in Fig. [THl 

The insulating behavior is accompanied by the change of the orbital 
ordering (the distribution of the Mo(4<i) electron densities). In the metallic 
state, it comes exclusively from the local trigonal distortions of the oxy- 
gen octahedra and represents the alternating a\g orbital densities in the 
background of degenerate orbitals (Fig. EDI). In the insulating state, the 
orbital ordering is determined not only by the local trigonal distortions, 
but also by the form of SE interactions between NN Mo sites, and tends 
to minimize the energy of these interactions [8j . Two typical examples for 
the FM and AFM (obtained after the flip of magnetic moments at the sites 
2 and 3) phases are shown in Fig. EH As expected for the FM spin order- 
ing [8j , the orbitals tend to order " antiferromagnetically" and form two 
Mo sublattices. Clearly, this orbital ordering breaks the Fm3d symmetry: 
if the sites belonging to the same sublattice can still be transformed to 
each other using the symmetry operations of the Fm3d group, the sites 
belonging to different sublattices - cannot. This will generally lead to the 
anisotropy of electronic properties, including the NN magnetic interactions. 
The AFM spin ordering coexists with the FM orbital ordering. It breaks 
the Fm3d symmetry in the spin sector, but not in the orbitals one. 

Results for NN magnetic interactions are shown in Fig.E2l as a function 
of U. We note the following. 
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Figure 19: Distribution of the Mo{Ad) states obtained in model Hartree-Fock calculations 
for C/=3.0 eV. 
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Figure 20: Orbital ordering in Nd2Mo207 for 17=1.5 eV. 




Figure 21: Orbital ordering in Y2M02O7 for C/=3.0 eV. Left panel: ferromagnetic phase. 
Two orbital sublattices are shown by black and gray colors. Right panel: antiferromag- 
netic phase. Two spin sublattices are show by black and gray colors. 

1. Jtt'j which are ferromagnetic for small [/, exhibit a sharp drop at the 
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Figure 22: Nearest-neighbor exchange interactions calculated in the ferromagnetic state. 
The orbital ordering realized for large U breaks the Fm3d symmetry, that leads to the 
inequality Ji2=Ji3i^Jia- Two such parameters, J12 and J14, are shown by dark and light 
symbols, respectively. 



point of transition into the insulating state. 

2. There is a significant difference between Nd/Gd and Y based com- 
pounds: in the case of Y, the exchange parameters are almost rigidly 
shifted towards negative values, so that the NN couphng become anti- 
ferromagnetic in the insulating phase, while it remains ferromagnetic 
in the case of Nd and Gd. 

The behavior can be easily understood by considering partial aig and 
e'g contributions to the NN exchange coupling, calculated after transfor- 
mation to the local coordinate frame at each site of the system (Fig. 12^ . 
Large FM e^-e^ interaction in the metallic regime is related with the dou- 
ble exchange (DE) mechanism, which is the measure of the kinetic energy 
for the itinerant f-spin electrons. As long as the system is metallic, the 
DE interactions are not sensitive to the value of [/, and the FM coupling 
will prevail. In the insulating state, the e'g electrons are localized at the 
atomic orbitals. This reduces the kinetic energy and suppresses the DE 
interactions, that explain the sharp drop of J^-t' (Fig. [221) • 

The main difference between Y and Nd/Gd based compounds is related 
with the aig-aig interaction. Since the |-spin aig band is fully occupied and 
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Figure 23: Contributions of aig and e'g orbitals (coming from the t2g manifold in the local 
coordinate frame) to the exchange constant J14. 

the J,-spin band is empty, this interaction is antiferromagnetic due to the 
superexchange mechanism. Since the SE couphng is proportional to the 
square of the aig bandwidth, this interaction will be the largest in the 
case of Y, that explains the AFM character of the total exchange coupling 
realized in this compound for large U. 

This example clearly shows at the importance of Coulomb U in the 
problem interatomic magnetic interactions of Mott insulators: it is abso- 
lutely indispensable in order to open the gap between occupied and unoc- 
cupied states, and suppress the FM DE interactions. Only in this case the 
total coupling between neighboring Mo sites can become antiferromagnetic, 
which is a necessary precondition for the SG behavior. On the contrary, 
the metallic state of Mo pyrochlores will always coexist with the ferro- 
magnetism. Therefore, it is not right to say that the energy gap is not 
a ground-state property and therefore does not need to be present in the 
Kohn-Sham quasi-particle spectrum of the SDFT. The present example 
shows exactly the opposite: the energy gap may determine not only the 
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value, but also the sign of interatomic magnetic interactions, which are the 
ground-state properties. 

6 Concluding Remarks 

In this review we tried to formulate some framework for the analysis of 
interatomic magnetic interactions in the transition-metal oxides and give 
some hints and ideas illustrating how this analysis should be generally 
done or at least started with on the basis of realistic electronic structure 
calculations. It is based on the fundamental magnetic force theorem, which 
allows to connect the problem of interatomic magnetic interactions in solids 
with the (generally auxiliary) single-particle spectrum obtained from the 
solution of Kohn-Sham equations in the spin-density-functional theory for 
the ground state. 

The straightforward application of this theorem is seriously hampered 
by the fact that in realistic calculations, SDFT is supplemented with some 
approximations, which are not always adequate for TMO. Nevertheless, in 
many cases the electronic structure calculations can provide the basic idea 
or at least some useful input information for the analysis of interatomic 
magnetic interactions in TMO, which (perhaps with some imaginations!) 
can be transformed to a realistic semi-quantitative model. This point of 
view was illustrated for a number of popular nowadays compounds, such 
as CMR manganites, double perovskites Sr2FeM06 (M= Mo, Re), and 
pyrochlores ^2Mo207 {A= Y, Gd, and Nd). 

It would be unfair to say that this strategy is universal and can be 
directly applied to all possible types of TMO. In some sense, we were 
lucky that the realistic information about the electronic structure of the 
considered compounds can be obtained already on the level of LSDA: 

• It gives a rough idea about the competition of double exchange and su- 
perexchange interactions in CMR manganites and double perovskites. 

• It allows to identify the electronic states, which are mainly responsible 
for the physics of Mo-pyrochlores and well separated from the rest of 
the spectrum, so that the next step connecting the LSDA calculations 
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with the degenerate Hubbard model becomes basically a matter of 
routine. 

Unfortunately, such a generalization is not always possible and there 
is a wide class of magnetic TMO, which still presents a very challeng- 
ing problem for the first-principles electronic structure calculations. This 
group includes many nickelates and cuprates. These systems are especially 
complicated because of two reasons: 

• For many of them, the standard LSDA calculations erroneously 
lead to the nonmagnetic ground state (the most famous example is 
La2Cu04 [77j). Therefore, it is difficult to make any conjectures bas- 
ing on the LSDA. 

• All these compounds are of the charge-transfer type. In this case, the 
semi-empirical LDA+f/ calculations may not be very reliable either. 

In order to proceed with the realistic description of electronic and magnetic 
properties of these materials, the development of new methods of electronic 
structure calculations, which would go beyond LSDA and get rid of several 
empirical features of the LDA+[/ approach is a very important task. 
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